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Universality of statistical properties of passive quantities advected by turbulent velocity fields at 
changing the passive forcing mechanism is discussed. In particular, we concentrate on the statistical 
properties of an hydrodynamic system with pressure. We present theoretical arguments and prelim- 
inary numerical results which show that the fluxes of passive vector field and of the velocity field 
have the same scaling behavior. By exploiting such a property, we propose a way to compute the 
' anomalous exponents of three dimensional turbulent velocity fields. Our findings are in agreement 

, within 5% with experimental values of the anomalous exponents. 
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Zero-mode dominance is of crucial importance because naturally explains the deviations from dimensional predictions 



I. INTRODUCTION 

Recent theoretical and numerical studies of passive quantities linearly advected by either stochastic j|,||,|l^,^,|8| or 
true turbulent Navier-Stokes velocity fields have focused on the important problem of the "zero-mode" dominance 
of the statistical properties of the passive field in the inertial range. 
The field may be an advected scalar 9{x,t), with the equation of motion 

f)0 

— +u-V9^K\7^e + fe, (1) 

or a vector, like a magnetic field B{x,t) satisfying M 
^ dB 

00 i — +uVB = B ■Vu + kV'^B +{b. (2) 

in ■ 
o 

for scaling properties and some universality with respect to the forcing mechanism, i.e. some independence of scaling 
y—( • properties from the large scale input of passive field. 

' In particular, for the special class of Kraichnan-like problems for scalar and vectors with Gaussian and white- in- 
time velocity field and external forcing |]l|,p|,p^P] , one can prove that the equal time (passive) correlation functions are 
dominated, in the inertial range, by the zero modes of the linear operator describing the advection by the Gaussian 
Jlj velocity field. Anomalous scaling comes from the non-trivial scaling properties of the null-space of the inertial operator. 
Universality comes from the fact that only prefactor of the zero modes feel the large scale boundary conditions, while 
the power law behavior is fixed by the inertial operator. As a consequence of zero-mode dominance in the inertial 
range we also have that all forcing-dependent contributions are sub-dominant and therefore some degree of universality 
. with respect to the forcing mechanism. 
■ ■ ■ '• Recently, numerical evidences that zero modes show up also in true passive models, advected by realistic non- 
Gaussian and non white-in-time velocity fields, have been presented for both passive scalars in the inverse cascade 
regime of a two dimensional fiow j|] and shell models for passive scalars [|j . Such findings give strong support to the 
idea that statistics of quantities linearly advected by true Navier-Stokes fields may show some universality properties 
with respect to the forcing mechanism too. In particular, the problem of a passive with a mean shear studied in |^ 
shows that zero-modes remain dominant with respect to the forcing contributions also in presence of a correlation 
between forcing and velocity advecting field, i.e. some degree of universality still holds. In this paper we want to 
investigate how far one can push the idea of universality with respect to the forcing mechanism for the particular case 
of the advection of a passive field with pressure advected by the true Navier-Stokes turbulent velocity field. Such 
a model was initially proposed in in order to understand the importance of non-local contributions, induced by 
the pressure term, in the zero-modes structure of the advecting operator. Here, we want to show and study another 
important striking feature of the model, i.e. its property to describe either a linear or non-linear evolution depending 
on the correlation between the passive forcing and the advecting velocity fields. Moreover, when the passive evolution 
is indeed non-linear, the passive field and the velocity fields coincide, i.e. the transported field and the advecting field 
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are looked in both their spatial and temporal evolution. In other words, we have a linear model which may become 
non-linear by changing the correlation properties between the external forcing and the advecting velocity field. 

In this paper we address two problems. First, we want to understand the degree of universality in the scaling 
properties of the linear problem at changing the external forcing. 

Second, by using the similarity between the non-linear Navier-Stokes eqs and the linear advecting-diffusive eqs for 
the passive field (see below) we propose a way to compute the anomalous exponents of homogeneous and isotropic 
turbulence which takes advantage of recent results in the theory of passive scalar advected by Gaussian, white in time 
velocity fields 

The paper is organized as follows. In section II we present the equations for the passive field with pressure, in section 
III we present some numerical experiments made on the equivalent shell model for the passive field and we discuss to 
which extent we may expect strong universality with respect to the forcing mechanism in this model. In particular 
we discuss why we expect, and verify numerically in the shell model case, that passive fluxes posses a much better 
degree of universality than passive structure functions. In section IV we used the, supposed, universality of passive 
flux in order to derive a functional constraint for the velocity statistics. Conclusions follow in section V. 

II. THE MODEL 

Let us consider the incompressible Navier-Stokes (NS) equations: 

dt + Uj dj U, ^-d,p + uAU, + f, (3) 

where notation follows the usual meaning and the forcing term fi is supposed to produce a stationary homogeneous and 
isotropic turbulence. We next consider a vector field W, divergence-less (i.e. diWi = 0) and satisfying a transport-like 
equation: 

dtW, + UjdjW, = -d,p^ + + 5, (4) 

where the "pressure" term is computed via the Poisson equation: 

d.UjdjW, = -Ap„ (5) 
Our aim is to discuss the statistical properties of W and their relation with the scaling properties of U . 

Let us introduce the anomalous exponents for U and W . We shall denote by C(p) the scaling exponents of the 
longitudinal structure function: 

where 5rU = {U{x r) — U{xy) ■ f and ( • ) stands for ensemble averaging. In the same way, we introduce (j{p) as 
the scaling exponents of longitudinal linear-fleld structure functions: 

where SrW = {W{x + r) — W{x)) ■ f. Finally, we consider the scaling exponents s{p) for the W flux defined as: 

{{5rUY{5rWfP) - r"(p) (6) 

Let us remark that equation (^) implies s(l) = which follows from the analogous of the "4/5" Kolmogorov equation 
for W Q, namely 

{5rU{5rWf) - Nr (7) 

where N is the mean rate of W dissipation. Our analysis is aimed at understanding the relationship, if any, among 
the anomalous exponents C(p), cr(p), s(p). We are now able to pose our problem in a quantitative way. Let us first 
consider the very simple case fi = gi. By subtracting equation (||) from equation (^ we obtain: 

dt(j)i + Ujdjcfii = -diTT + L'A(j>i (8) 
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where (p = W — U and tt = {p — p^). By equation (||) it immediately follow that the space average i?^ of 0^ = 
J dx(j)i{x)(j)i{x) satisfies the equation: 

dtE^ = (9) 

where is the mean rate of dissipation of E^. Thus, for long enough time, the field c/fi goes to zero and Wi — Ui, 
identically. 

In this paper we want to understand to which extent the statistical properties of are universal with respect to 
the forcing mechanism. If a strongly universality holds, then we should have independently of the forcing mechanism 

cT{p)^ap), s(p) = c(3p) (10) 

even when fi and gi are uncorrelated or weakly correlated fields. Let us notice that in the previous equality we have 
assumed that velocity flux possesses the same scaling properties of the velocity field, as always verified in all numerical 
and experimental data. We do not know any rigorous argument against or in favor of universality for the scaling 
properties of (^). One may argue that if the forcing mechanisms fi and gi are weakly correlated -or independent- 
and confined only to large scales than the statistical properties of W may not be strongly influenced by the forcing 
itself. In the latter case the universality of the scaling exponents of W should be achieved by the same zero-modes 
mechanisms previously discussed for the passive scalar with independent forcing. Whether this supposed universality 
can be pushed until the very extreme case of fully correlated systems, implying the equalities ( p^ ) is a matter of 
discussion. There is clearly a physical relevant question here: we need to understand the importance of external 
forcing mechanism in the statistical properties of the advected field. This question, of course, arises not only for the 
passive vector with pressure but also in any linear advection problem, i.e. passive scalars or passive magnetic fields. 
In the case of the passive vector with pressure the question assumes also another interest because of the possibility 
to push the advected field to follow the velocity field for some particular external forcing. 

We consider an important point to investigate this problem in details and in the next section, we present some 
numerical results showing that universality holds for fluxes, i.e. s{p) — C{3p). Passive structure functions seems to 
be less universal than fluxes although strong boundary effects (both ultraviolet and infrared) do not allow a precise 
statement. 



III. NUMERICAL RESULTS 

Although a direct numerical tests of the statistical properties (0) is possible for the NS equations, we will limit in 
this paper to a direct numerical investigation in the framework of shell models for three dimensional turbulence. 
Shell models for turbulent energy cascade have proved to share many statistical properties with both turbulent three 
dimensional velocity fields and with passive linearly advected quantities Let us introduces a set of 

wavenumber fc„ = A^/cq with n = 0, . . . , N and the inter-shell ratio fixed to A 2. The shell- velocity variables u„(t) 
must be understood as the velocity fluctuation over a distances r„ = k~^. We also introduce the linear- field shell 
variables Wn{t) defined on the same sets of discrete wavenumber. It is possible to write down a set of coupled ODEs 
for the time evolution of Un{t) and Wn{t) which mimics the velocity and the passive turbulent evolution (see below). 

We know that zero modes are at works in shell models for linearly advected quantities, exactly like in the true linear 
hydrodynamical problems. Indeed, it is possible to prove analytically that shell models for passive scalar advected 
by Gaussian and white-in-time velocity (shell) fields have intermittent corrections dominated by the null space of 
the linear finite-dimensional advecting operator Recently, it has also been shown that shell models for passive 

scalar advection, i.e. passive shells advected by shell fields arising from a shell model for the velocity, have zero modes 
dominance for the scaling properties in the inertial shells, exactly as it is shown for the true passive scalar [^. 
We have therefore, exactly the analytical/phenomenological framework useful to check to which extent the scaling 
properties of the linear hydrodinamical problem are forcing independent. 

In the shell model framework the very meaning of pressure is absent. The equivalent of our linear- hydrodynamical 
model (Q) will become a shell field Wn{t), linearly advected by the shell model velocity fields Un{t), conserving the 
energy l^np, and such that when /„ = gn we have, for time large enough, Wn{t) = Un{t). 

It is possible to write down such a coupled set of ODEs for all shell models. Here we consider the case of Sabra 
model |l|]. We obtain: 

(-^ + '^kl)Un = i(fc„U*+iM„+2 + bkn-lUn+lUn-1 + (1 + 6)A:„_2U„_2U„_i) + /„ (11) 
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+ vk'^i)Wn = i(A:„M*+ilO„+2 + &fc„-lW„+lM*_i + kn-2Wn-2Un-l + 6fcn-2W„_iU„_2) + ffn (12) 

where the non hncar term of it„ evolution and the hnear advection part of Wn evolution have the only free parameter 
b. Note that if we put ui„ = m„ in equation ( p^ ) , then equation (^) becomes equivalent to equation ( pi] ) . It is known 
that in order to have a realistic intermittent behavior for m„ one has to chose — 1 < 6 < jlj]. 

Equations (^l|) and ( p^ ) preserve velocity energy, = I'^nP' passive energy, E^j = I'^^nP, in the limit of 
zero viscosity and zero forcing. 

It is easy to show that if /„ = gn in the above equations than for long enough time Wn{t) = u„(t). 

The universality issue, as discussed in the previous section, consists now in studying the scaling properties of w„ 
at changing its forcing mechanism. 

Let us define the scaling exponents for the velocity flux, 

n„ = 5[(fc„7i„+2Mjl+iU*J + (1 + b)kn-l{Un+lU*^U*_i)] 

and passive flux, 

Qn = 3[(/c„w„+2W*+iU'*) + fc„_i(w;„+iu*w*_i) + fc„_i6(w„+iw*u*_i)] 

as: 

where the equivalent of 4/5 law for the two shell models gives C(3) = s(l) = 1, The structure functions are 

defined as follows: 

SZin) = {\u,,n ^ k-<^P^ TP{n) = {\w^n ^ k'^^^^ (14) 

In the following we present some numerical tests done with = 25 shells Uu — ~ 5 ■ 10^^ and different kind of 
forcing mechanisms. 

Forcing have been chosen such as to go from a weakly correlated situation where gn (t) has some large scale dependency 
from the Un(t) dynamics to a fully uncorrelated case with gnit) given by a random process. 

In Fig. 1 we present the typical scaling laws one obtains for fluxes of both fields, for different large scale passive- 
forcing. In all simulations we have always taken the same velocity forcing concentrated on the largest shell and 
constant: = (1 + i)C„(5„,i, with C„ — 0.01. We have compared statistical properties of the passive fields using 
three different choices for the passive forcing. In case (A) we had a time-dependent forcing such as to impose a 
constant passive energy input on the first shell: 

A: gn{t)^6n,i{'^ + i)/wl{t). 

In case (B) we fixed the first passive shell to have the same amplitude Itfil = const, but leaving its phase to evolve 
according to its own dynamics. 

B : gn{t) -> \wi{t)\ = 1. Vi. 

In case (C) we took a forcing concentrated on the first shell, with constant amplitude, |(?i(i)| — Gi but random 
independent phases 

C: gn{t) ^ Sn,iGie''^'^ 

where {e{t)6{t')) oc S{t - t'). 

Let us notice that forcing of cases (A) and (B) have some (weak) correlation with the adverting velocity field, while 
case (C) is independent of u„. 

As one can see all flux curves superpose perfectly in the inertial range. Let us notice that the extremely small errors 
on the scaling of two fluxes allows us to support the statement n„ ^ (5„ with very high accuracy independently on the 
forcing mechanism. The scaling of passive structure functions suffers of larger error bars, due to a less smooth matching 
between inertial and infra-red properties [ p6[ . The qualitative trend is, anyhow, toward a more intermittent behavior 
of the passive fields with respect to the velocity field (see Fig. 2). In the latter case we would have an indication that 
passive structure functions are much more sensible to the boundary conditions than the fluxes moments, implying 
also a weaker degree of universality with respect to the forcing mechanism. 
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In Table 1 we quantify our finding in all three forcing cases (A-C) showing the best fits and their errors for all scaling 
exponents of both fluxes and structure functions. 

In Fig. 3, we show two typical time evolution for the total energies of velocity, Eu{t), and passive, Eu,{t) for the case 
with two different -but correlated- forcings, /„ ^ gn- Although some correlation between Eu{t) and E^{t) is observed, 
we are very far from the trivial exact correlated case one would have obtained choosing /„ = g-a- The non trivial 
correlations between the two energy can be seen in the fig. 4 where we plot, for a typical trajectory, E^ versus E^. 
In order to test the robustness of previous results we repeated the numerical experiments with a different value of 



the free parameter h in the sabra shell model equations (11). It is known that at changing h the intermittency of the 
velocity field changes. 

In Fig. 5 we plot the C(3p) and s{p) curves for both values of 6 = —0.4 and b — —0.6. The agreement for each h values 
between the two fluxes is again perfect despite the fact that at varying b we have different degrees of intermittency. 

One could question why there is such clear evidence that s{p) — C(3p), i.e. the scaling properties of two fluxes are 
identical, while passive structure functions seem to be more intermittent than the corresponding velocity structure 
functions. One possible explanation goes as follows. By exploiting the equation of motion one may always derives 
homogeneous constraints for moments of quantities like En — UnWn'Wn" in the inertial range, i.e. where forcing is 
not directly acting. In particular by writing the stationary condition for quantities like (Fni^n, • • • En^ -^WkWk') one 
obtains an homogeneous constraints involving only Fn observable: 

J^iFn.Fn. ■ ■ ■ Fn,_,w,w,,) = d^^^^r^l^^^,(Frn.Frn. ■ ■ ■ Em,) = (15) 

where the operator OjJ^^^J^.^"^^ ^ ^, is given from the equation of motion and it is independent of the correlation 
between m„ and w„ and of the chosen forcings gn and /„. On the other hands no homogeneous closed constraints 
can ever be found for moments of quantities involving only simultaneous correlations of passive fields, w„. In order 
to have homogeneous expressions for obervable made of only passive fields one has to give (or to guess) an explicit 



form for the correlation between w„ and u„. Thus one may think to obtain an equation similar to (15) but now 
with the inertial operator explicitly dependent on the given correlation between the two fields. In the latter case 
the dependency of the passive field on the statistical properties of the velocity field and of the correlation with the 
external forcing may lead to non-universal scaling properties. 

On the other hand, one can also argue that the observed lack of universality is not Reynolds independent and that 
at Reynolds large enough some strong independence from the forcing mechanism (zero-modes dominance) would 
be recovered also for passive structure functions. The latter scenario is what happens in true passive scalars with a 
correlation between forcing and velocity field (see the case of a passive scalar with shear |^) where the existence of sub- 
leading non-homogeneous terms induced by the forcing mechanism may spoil the scaling behavior of the zero-modes 
at finite Reynolds numbers. 



IV. AN APPROXIMATE COMPUTATION OF THE SCALING EXPONENTS 

The results so far discussed, make us rather confident that the two fluxes of equation (||) and of Navier Stokes eqs 

have the same statistical fluctuations. 
We want now to understand whether it is possible to use the above results to obtain useful informations on the 
scaling exponents C(?')- 

The main idea, discussed in this section, is the following. We assume that the velocity field can be described by 
the (unknown) multifractal probability distribution P{5rU). We want to compute the probability distribution of the 
passive vector, P{5rW). The equation (^) will induce a functional relation between the two probability distributions. 
By requiring that the scaling properties of W flux are the same of U flux, we introduce an (infinite) set of equations 
for the probability P{5rU) whose solutions will fix its scaling. In order to simplify the above procedure, the approach 
we want to follow is based on a suitable set of approximations and assumptions whose validity we are not able to 
prove rigorously, only a posteriori we will be able to judge the goodness of our calculation. 

Let us suppose that the field U in equation (^) is characterized by a multifractal spectrum D{h). As it is well known, 
in such a case the anomalous exponents of U are given by the expression C(p) = iiiih{ph + 3 — D(h)) and the probability 
distribution of 6rU is given by [T^ : 

P{SrU) (X / dn{h)r^-^^'''> cxp 



(16) 
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where Uq represents the variance of the large scale velocity field, which is supposed to be Gaussian. Equation ( |l^ ) 
tells us that for each value of r we may consider SrU as the superposition of Gaussian field with variance Uor^'^ and 
probability r^-^W associated to each value of h. We can then consider to solve equation (Q) in the limit where U is 
a "weighted" superposition of Gaussian random fields (as expressed by (Of 



The equivalence between the statistical properties of the Navier-Stokes field and of the linear field implies that both 
fiuxes have the same scaling: 

{\{SrU)iSrW)^\n {\SrU\'n (17) 

We want to exploit the above identity in order to derive a constraint for the D{h) spectrum. 
We can rewrite the term {{SrU{SrW)'^)''') as follows: 

{(SrUiSrWff) ^ I dhdiSrUy~^^''^P{6rU,h)i6rU)P{i6rWfP\6rU,h} (18) 



where with {{6rW)'^P\6rU, h) we mean the average of the linear field conditioned to the advecting velocity field and 
P{6rU, h) = exp — 260)^*^ ■ ^l"^' ( P^ '^ot nsed any approximation. The computation of {{5rW)'^P\5rU, h) 

is the most difficult one and we shall introduce several approximations in order to make it feasible. In particular we 
want to compute {{5rWYP\5rU,h) using a surrogate (5-correlated and Gaussian velocity field. In other words will 
make the first approximation that in order to compute the conditional average we may assume that the multifractal 
velocity field, as defined by equation (|l6|), is the "weighted" superposition of independent Gaussian random fields 6 
-correlated in time. 



In order to compute left hand side of equation (18), we solve equation (^ for fixed exponent h of the random field U 
and then average the results over h with probability r^-^W . In doing such an average we notice that in order to mimic 
the advection by a true Navier-Stokes field (not J-correlated in time) with scaling 5rU ~ we need a J-correlated 
surrogate with scaling SrUg ^ r~5~, because of dimensions carried by the delta-functions. As a consequence we are 
going to be interested only in exponents for the surrogate Us field in the range H = = [1/2, 1], which correspond 
to the range h : [0, 1] for the true turbulent field. Let us notice that a K41 field has a correspondent (5-correlated field 
scaling as 5rUs r^^^. 

We now want to exploit our (5-correlated ansatz by writing: 

{{SrWfP\SrU, h) - rf(i-'')+''^'' (19) 

where we have introduced the anomaly p2p(2H) of the linear field W advected by a (5-correlated velocity field Srlls ~ 
r^, namely: 

{{SrWfn)H rf(2-2H)+p.,(2H) (20) 

In equation ( pO| ) we have taken into account the dimensional consistency relation 2H = 1 + h. Notice that for the 
(5-correlated case one can prove that in the isotropic sector we have p2{'2H) = 0, i.e. the second order correlation 
function has not anomalous scaling. 
Our definition gives: 

{5rU{SrW)^)H - r"r^-" = r (21) 

i.e, our ansatz on the function {{5rWY'P\5rU , h) implies that we are averaging over all possible singularity H with the 
constraint {5rU{5rWY) ^ r. 

We can now use the relation of statistical identity between fiuxes, C(3p) = ^(p)? i^i order to obtain an equation for 
D{h). We have: 



which gives: 



d^(/i)r^'''+P(^"''^~''"''(^+'')+^"-°(''' - /^^P' (22) 



inf [p - p2p(l + /i) + 3 - D{h)] = inf [dph + 3 - D{h)] (23) 

h h 
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If we are able to compute the anomalous correction p2p (of the (5-correlated problem), equation (^3|) becomes a 
functional equation for D(h) whose solution gives us the anomalous exponents C(p) within the set of approximation 
discussed above. Let us notice that ( p^ ) is consistent with the known result r(l) = s(l) = 0. Let us also remark that 
as soon as we introduce any smoothing in the time dependency of U the equation ( p3| ) may be no longer valid and 
we should reconsider the averaging procedure in a more suitable way. 

The computation of the anomalous exponents p2p is a feasible but difficult task for equation because of the 
condition diWi = 0. In principle, the computation can be done perturbatively following the analog of the passive 
scalar case ||^. Only the exponents in different anisotropic sectors for the second order correlation function have 
been, up to now, computed. However, in order to understand the quality of the approximations so far introduced to 
derive equation (p3|), at least for three dimensional isotropic and homogeneous turbulence, we can use the numerical 
values oi pi{\ + h) and pq{\-\- h) as recently computed numerically for the passive scalar model [ p^ , hoping that the 
introduction of the pressure term does not change too much the scaling exponents pO| . 

In order to find a solution of (^^, we simply assume that D{h) can be parameterized by using the expression 
D(h) = do(l — a; + xl'n{x)) where x = {h — h^) / (dolog{l3~^)) which corresponds to a log-poison distribution for 
the multifractal model. We use a log-poison formula because we know that it parametrize particularly well the 
experimental data ||2l| , |2^ . Because do is fixed by the condition C(3)_= 1, we are left with the problem to find /? and 
ho by using equation j23|) for p = 2 and p = 3. Solving equation ( [2^ ) for /? and ho gives /3 = 0.781 and Hq = 0.14. 
In fig. 6 we show the estimated values of ({p) for p = 1, ..12 against the experimental findings for homogeneous and 
isotropic turbulence. We note that the estimate based on ( |2^ ) is rather accurate with an error not greater than 5%. 
We argue that such a small error shows that our approach, within the limitation and the approximations previously 
discussed, looks promising as an useful tool to compute the anomalous exponents. 

In order to test the sensitivity of (23) to the values of p4{l + h) and pQ{l + h), it is interesting to use the anomalous 
exponents pi and pe given by Kraichnan' formula for D = 3. In this case the solution of ( p3| ) gives /3 = 0.813 and 
ho — 0.21 and the corresponding values of ({p) are also plotted in fig. 6. Although there are not big differences 
between exponents obtained using the Kraichnan formula and those obtained by using the correct numerical results, 
the agreement with the experimental data is definitely better in the latter case. 



V. CONCLUSIONS 



In this paper we discussed several properties characterizing the statistical behavior of a divergence-less vector field 
passively advected by an homogeneous and isotropic turbulent field. We can summarize our findings in the following 
way. (i) We present some theoretical arguments which supports the statement that the flux of passive vector field 
should have the same statistical properties of the flux of an homogeneous and isotropic turbulent field, (ii) We 
generalize the concept of divergence-less vector field to shell models and by using detailed numerical simulations we 
provide evidence that the anomalous scaling exponents of the passive vector fiux and of the non linear shell model flux 
are the same with very high accuracy, (iii) We propose a self consistent approach to compute the anomalous exponents 
in homogeneous and isotropic turbulence by using the properties previously discussed. In particular, we have been 
able to define a functional equation for D{h) within a suitable set of approximations, (iv) By assuming that the 
pressure term does not change dramatically the numerical values of the anomalous exponents for the divergence-less 
vector field advected by a Gaussian isotropic (5-correlated random field, we have been able to find an approximate 
solution of the functional equation for D{h) which compares rather well to known experimental data. 

Before closing we want to discuss a generalization of the model (^ which we consider interesting for further studies. 
Let us consider the following set of equations: 

dtU, + U.VjU, + XW.VjU, = -V^pu + i^AU, + f, (24) 



dtW^ + UjWjW, + XWjVjW, = -V,p^ + vAW, -f g, (25) 

where we assumed that dilJi ~ diWi = 0. It is interesting to notice that the vector field Z\ = U + XW satisfies the 
Navier Stokes equations and that equations (p^ and (|2^) corresponds to (|^) and (^ for A = 0. 

Finally we want to remark that equation (0)^ is an useful tool to understand the role played by coherent structure 
on the anomalous scaling in both two dimensional and three dimensional turbulence. All our findings suggest that 
a systematic study of equation (^ looks extremely promising in order to derive a new approach in understanding 
intermittency and estimating anomalous scaling in Navier-Stokes equations. 

We acknowledge useful discussion with I. Arad, A. Celani and M. Vergassola. We thank M. Vergassola and A. Mazzino 
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and by the INFM (Iniziativa di Calcolo Parallelo). 
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Jr* 


IV7J, n \ / 


fe-f/^r^/^ra;) (A) 


kZ^^'-^Tl''-^in) (B) 

Q V") \ ) 


k^^''-^Tl''-^in) (CI 

^ Q v / \ ) 


TK(n) (A) 


TE.in) (B) 


TPAn) (C) 


2 


0.712 (3) 


0.711 (3) 


0.711 (2) 


0.710 (2) 


0.67 (3) 


0.67 (3) 


0.66 (3) 


4 


1.263 (6) 


1.264 (7) 


1.264 (6) 


1.266 (3) 


1.20 (6) 


1.28 (4) 


1.13 (7) 


6 


1.745 (8) 


1.741 (7) 


1.74 (1) 


1.745 (5) 


1.6 (1) 


1.67 (8) 


1.4 (1) 


8 


2.18 (2) 


2.18 (2) 


2.19 (2) 


2.18 (2) 


2.0 (2) 


2.1 (1) 


1.6 (2) 


10 


2.60 (2) 


2.58 (3) 


2.61 (2) 


2.57 (4) 


2.3 (2) 


2.6 (2) 


1.9 (3) 



TABLE I. Scaling exponents of both flux and structure function of the passive field and of flux of velocity fleld at changing 
the order of the moment p = 2,4,6,8,10 and for different large scale forcing cases (A-C). Notice the high precision in the 
agreement between the passive flux at changing the large-scale forcing (columns 1-4). Passive structure functions (columns 
5-7) are less accurate due to the presence of a strong bottleneck at small scales. Errors are given by the numbers in brackets 
and refers to the last digit. Errors are estimated from the fluctuations of the local logarithmic derivatives in the shell range 
n = 3, . . . , 15. 
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FIG. 1. Log-log plot of velocity and passive fluxes and structure functions for second order moment, p = 2 versus the scale 
kn- Curves represent from above: (□) velocity flux S^(n). Passive flux, Tqin), for forcing case (A) (x); forcing case (B) (+); 
forcing case (C) (*). Fluxes are always multiplied by the normalising factor fc^^. Curves have been shifted along the y-axis for 
the sake of clarity 
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log2{kn) 

FIG. 2. Log-log plot of second order velocity flux, p = 2, and sixth order passive structure functions, versus the scale fc„. 
Curves represent from above: (□) velocity flux Suin); (x) passive structure function T^(n) forcing case (A); (+) passive 
structure function T^{n) forcing case (B); (*) passive structure function T^(n) forcing case (C). Notice the strong ultraviolet 
and infrared effects in the passive structure functions. Notice also that for the fully independent forcing, case (C), the passive 
structure functions shows a larger intermittent slope at small scales. 
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FIG. 3. Typical time evolution of the total velocity energy Eu{t) (solid) and total passive energy E,j,{t) (dashed) (case A). 
Curves have been shifted along the y-axis for the sake of clarity 
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FIG. 5. Scaling exponents of velocity flux and passive flux for two different values of the parameter b in the "Sabra" 
shell model equations (|ll|). □ corresponds to the velocity flux exponents, (^(p/3), and the dashed curves to the passive flux 
exponents, (j{p). Above: the case with b = —0.4, below the case with h = —0.6. Notice that despite the remarkable change in 
the intermittent properties the fluxes follow eachother perfectly. 
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FIG. 6. Comparison between the scaling exponents computed by eqn. (ESI) using: numerical results for the passive scalar 
[17] (solid curve); the Kraichan formula for the passive scalar [1] (dashed curve). Experimental data (+) are given with their 
error bars. 
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